Stratum dynamic correction method of guidance while drilling

专利摘要:
The present disclosure relates to a stratum dynamic correction method of guidance while drilling. The method comprises the steps of obtaining an expert model and layered data points of stratum interpretation; adding control points according to the real-time elevation data of the expert model; according to the layered data points and the control points, calculating elevations of grids influenced by all the control points, and correcting the expert model according to the elevations of the grids influenced by all the control points. According to the method, the models are selectively updated according to the grid points of the nearby influence areas where the control points are located, elevation values of all the grid models do not need to be reconstructed, the calculated amount is reduced, the updating speed of the models is greatly increased, and the real-time performance and precision of expert models are improved.
公开号:NL2028480A
申请号:NL2028480
申请日:2021-06-18
公开日:2021-09-16
发明作者:Wang Zhonghao；Liu Shaohua
申请人:Univ Yangtze；
IPC主号:

专利说明:

[01] [01] The present disclosure relates to the field of geological exploration, relates to a correction method of a geological model, in particular to a stratum dynamic correction method of guidance while drilling.BACKGROUND ART
[02] [02] In the interpretation process of real-time well logging while drilling, interpretation experts can obtain the expert model of stratum interpretation through the interpretation of the relationship between a borehole trajectory and a stratum, and make fine adjustments to the current stratum interface through this model. At this time, it is necessary to locally correct the initial structural model established before, so that the modeled stratum model is consistent with the actual stratum, providing a more accurate stratum interface for subsequent geological guidance and providing technical support for guidance.
[03] [03] Elevation refers to the distance from a certain point to absolute datum in the direction of the plumb line, which is referred to as absolute elevation and is simplify referred to as elevation. The distance from a certain point to the datum of an assumed level is referred to as the assumed elevation.
[04] [04] Horizon refers to a specific position in the stratigraphic succession. There are many stratum horizons, such as lithologic horizons with special lithologic characteristics, fossil horizons with special fossils, chronological horizons with specific era, seismic horizons, electrical logging horizons, etc. Therefore, the stratum horizon can be the boundary line of stratum units, and can also be a marker layer belonging to a specific era.SUMMARY
[05] [05] Aiming at the technical problems existing in the prior art, the present disclosure provides a method for drilling stratum comprising applying a dynamic correction guidance while drilling, so as to improve the precision and real-time performance of an expert model.
[06] [06] According to the present disclosure, the technical scheme for solving the technical problems is as follows. A method for drilling stratum comprising applying a stratum dynamic correction of guidance while drilling comprises the steps of: S1, obtaining an expert model and layered data points of stratum interpretation; S2, adding control points according to the real-time elevation data of the expert model, S3, according to the layered data points and the control points, calculating elevations of grids influenced by all the control points; and S4, correcting the expert model according to the elevations of the grids influenced by all the control points.
[07] [07] In some embodiments of the present disclosure, the step of calculating the elevation in S3 is as follows: calculating the distance from the current grid point to the layered data point and the control point in the area influenced by the control point; according to the distance between the current grid point, the layered data point and the control point, determining the elevation weights of the layered data point and the control point; according to the elevation of the layered data point and the control point and their respective elevation weights, obtaining the elevation of the current grid point.
[08] [08] Common distances include at least one of Euclidean distance, Manhattan distance, Chebyshev distance, Hamming distance and Minkowski distance. Further, the distance is Euclidean distance.
[09] [09] In some embodiments of the present disclosure, the distance between the grid point and the layered data point and the control point calculates the elevation weight, such as Euclidean distance. The elevation weight is inversely proportional to the distance from the current grid point to the layered data point and the control point.
[10] [10] In some of the above embodiments, the calculation method of the elevation weight is as follows: bos | dy jl a;
[11] [11] :
[12] [12] where pi is the weight of the i-th control point to the current grid point gridirow jet, di 18 the distance from the i-th control point to the current grid point gridiow jcot, and dj is the distance from the j-th control point to the current grid point gridiow jeol.
[13] [13] In some embodiments of the determined elevation weight calculation method described above, the calculation method of the elevation value of the current grid point is as follows: Spore goad > > Pry
[14] [14] Fd
[15] [15] where irow and jcol are the row number and column number of the previous grid, Ziow jcot ís the elevation of the current grid point gridicow jeot, pris the corresponding elevation weight of the control point or the layered data point, and z; is the elevation of the control point or the layered data point.
[16] [16] In some embodiments of the present disclosure, the step of calculating the elevation further comprises stopping the calculation when the elevation change of the current grid point is less than a threshold value. Preferably, the threshold value is not more than 0.01m.
[17] [17] In the above embodiments, the real-time elevation data all include the elevation change determined by actually measured borehole data, profile data, plane data and result data recognized by experts.
[18] [18] The method has the beneficial effects that the models are selectively updated according to the grid points of the nearby influence areas where the control points are located, elevation values of all the grid models do not need to be reconstructed, and the updating speed of the models is greatly increased.BRIEFT DESCRIPTION OF THE DRAWINGS
[19] [19] FIG. 1 is a basic flow chart of a stratum dynamic correction method of a geological model according to some embodiments of the present disclosure;
[20] [20] FIG. 2 is a schematic diagram of a three-dimensional model of a stratum according to some embodiments of the present disclosure;
[21] [21] FIG. 3 is a schematic diagram of an adjusted three-dimensional model of a stratum according to some embodiments of the present disclosure.DETAILED DESCRIPTION OF THE EMBODIMENTS
[22] [22] The principles and features of the present disclosure will be described hereinafter with reference to the accompanying drawings. The examples are only used to explain the present disclosure, rather than limit the scope of the present disclosure.
[23] [23] The present disclosure relates to a stratum dynamic correction method of guidance while drilling, wherein: the method comprises the steps of: obtaining an expert model and layered data points of stratum interpretation; adding control points according to the real-time elevation data of the expert model; according to the layered data points and the control points, calculating elevations of grids influenced by all the control points; and correcting the expert model according to the elevations of the grids influenced by all the control points.
[24] [24] In some embodiments of the present disclosure, the step of calculating the elevation is as follows: calculating the distance from the current grid point to the layered data point and the control point in the area influenced by the control point; according to the distance between the current grid point, the layered data point and the control point, determining the elevation weights of the layered data point and the control point; according to the elevation of the layered data point and the control point and their respective elevation weights, obtaining the elevation of the current grid point.
[25] [25] Common distances include Euclidean distance, Manhattan distance, Chebyshev distance, Hamming distance and Minkowski distance. In order to improve the model precision, one or more of the distances can be selected for combined calculation. Preferably, the distance from the current grid point to the layered data point and the control point is Euclidean distance.
[26] [26] In some embodiments of the present disclosure, the elevation weight is inversely proportional to the distance from the current grid point to the layered data point and the control point. Specifically, the calculation method of the elevation weight is as follows: oS]
[27] [27] © formula (1)
[28] [28] where pi is the weight of the i-th control point to the current grid point gridirowjcot, di is the distance from the i-th control point to the current grid point 5 gridiowjcot, and dj is the distance from the j-th control point to the current grid point gridisow col.
[29] [29] In some embodiments of the determined elevation weight calculation method described above, the calculation method of the elevation value of the current grid point is as follows: 3 i. Sirew geel > > Pry 0 [30] Fed formula (2)
[31] [31] where irow and jcol are the row number and column number of the previous grid, Zirow jcot Is the elevation of the current grid point gridicow,jcot, Px is the corresponding elevation weight of the control point or the layered data point, and z; is the elevation of the control point or the layered data point.
[32] [32] In some embodiments of the present disclosure, a stratum dynamic correction method of guidance while drilling comprises the following specific steps:
[33] [33] (1) adding three-dimensional control points in the target horizon;
[34] [34] (2) calculating the distance between the current grid point gridiow jcot in the area influenced by the control point and the layered data point of the target layer;
[35] [35] (3) calculating the distance between the current grid point gridiow jc: and the control point;
[36] [36] (4) calculating the weight according to the above formula 1;
[37] [37] (5) calculating the elevation value of the current grid point gridiow jest according to the above formula 2;
[38] [38] (0) circulating all grids in the area influenced by control points, and repeating the above steps 2-5 to complete the adjustment of the model.
[39] [39] In some embodiments of the present disclosure, the step of calculating the elevation further comprises stopping the calculation when the elevation change of the current grid point is less than a threshold value. Preferably, the threshold value is not more than 0.01m.
[40] [40] Some embodiments of the present disclosure will be explained with reference to the attached drawings. The red dots P1 to P8 in FIG. 2 are layered data control points. interpretation experts can obtain an expert model of stratum interpretation by interpreting the relationship between a borehole trajectory and a stratum, and add a new control point P9 from elevation data. When the point pointed by the arrow in FIG. 3 is the layered data point which is the added stratum control point, the model needs to be adjusted, but it does not need to be re-modeled, and only the elevation of the grid points near P9 needs to be recalculated.
[41] [41] According to the coordinates of point P9, the row and column number where it is located are calculated, diffusing around the grid, and increasing outward step by step at a time. The elevation of the current grid point gridiowjcot is calculated. First, the distance between the grid point gridiow jot and points P1 to P9 is calculated, the weight pi is calculated by formula (1), and then the elevation of the current grid point gridirow jcot is calculated by formula (2). After outward diffusion, the change (absolute value of the difference) between the calculated elevation of the grid point and the original elevation 1s less than a threshold value, and the threshold value is not more than 0.01m. The outward diffusion is stopped. At this time, the model adjustment of the area influenced by the new control point P9 ends, and the whole model completes updating. The updated model is shown in FIG. 3.
[42] [42] In the above embodiments, the real-time elevation data all include the elevation change determined by actually measured borehole data, profile data, plane data and result data recognized by experts.
[43] [43] The above is only a preferred embodiment of the present disclosure, and is not intended to limit the present disclosure. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present disclosure shall be included in the scope of protection of the present disclosure.

权利要求:
Claims (10)
[1]
A method for drilling stratum comprising applying a stratum dynamic correction of the guidance during drilling, wherein: the method consists of the following steps: obtaining an expert model and layered data points of stratum interpretation; adding checkpoints according to the real-time evaluation data of the expert model; calculating heights of grids affected by all control points according to the layered data points and the control points; and correcting the expert model based on the heights of the grids affected by all control points; and adjusting the borehole trajectory of the drilling equipment according to the corrected model.
[2]
The method of claim 1, wherein the step of calculating the height is as follows: calculating the distance from the current grid point to the layered data point and the checkpoint in the zone affected by the checkpoint; determining the height weights of the layered data point and the control point according to the distance between the current grid point, the layered data point and the control point; obtaining the height of the current grid point according to the height of the layered data point and the control point and their respective height weights.
[3]
The method of claim 2, wherein the distance is at least one of Euclidean distance, Manhattan distance, Chebyshev distance, Hamming distance and Minkowski distance.
[4]
The method of claim 3, wherein the distance is Euclidean distance.
[5]
The method of claim 2 wherein the height weight is inversely proportional to the distance from the current grid point to the layered data point and the control point.
[6]
The method of claim 4, wherein the height weight calculation method is as follows: the Td where pi is the weight from the ith control point to the current grid point, rastefrowjcot, ie the distance from the ith control point to the current grid point, rasteTiowjcot, and dj is the distance from the jth control point to the current grid point, rasterirow Jeol.
[7]
The method of claim 6, wherein the calculation method of the height value of the current grid point is as follows: UR, Sirofoal = > PES; js where irow and jcol are the row number and column number of the previous grid, Ziow jor is the height of the current grid point, rasteriowjco, px is the corresponding height weight of the control point or layered data point, and z; is the height of the control point or layered data point.
[8]
The method of claim 2, wherein the step of calculating the height further comprises stopping the calculation when the height change of the current grid point is less than a threshold value.
[9]
The method of claim 1 or 2, wherein the threshold value is not greater than 0.01 m.
[10]
The method of claim 1 or 2, wherein the real-time evaluation data all include the elevation change determined by actually measured borehole data, profile data, face data, and result data recognized by those skilled in the art.

类似技术:

---

公开号 | 公开日 | 专利标题

---

CN105807316B|2018-04-03|Ground observation microseism velocity model corrections method based on amplitude superposition

---

WO2009062935A1|2009-05-22|Forming a geological model

---

CN105317375B|2017-10-17|Horizontal well is inducted into Target process and device

---

CN105093319B|2017-11-28|Ground micro-seismic static correcting method based on 3D seismic data

---

CN105631753B|2020-03-10|Horizontal well oil reservoir profile modeling method based on stratum attitude

---

CN104090301B|2017-08-25|A kind of method for asking for three-D high frequency static correction value

---

CN104730581A|2015-06-24|Method for locating microseism event point

---

CN102879820A|2013-01-16|Three-dimensional surface layer model construction method based on triangular grid

---

CN106295873B|2017-10-03|A kind of Grey Markov Chain trajectory predictions method towards coal cutter memorized cutting

---

CN111060135A|2020-04-24|Map correction method and system based on local map

---

RU2471152C1|2012-12-27|Method of aircraft navigation

---

NL2028480A|2021-09-16|Stratum dynamic correction method of guidance while drilling

---

CN108333627B|2019-01-08|Igneous rock area is broken recognition methods and the device of the true and false

---

CN107843928A|2018-03-27|A kind of shot point method for correcting error based on grid-search algorithms

---

US6016461A|2000-01-18|Method for producing maps of the risk in positioning a well in an environment

---

CN105911592B|2017-11-10|A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint

---

US6446006B1|2002-09-03|Method of producing impact point position probability maps for a well

---

KR20200058258A|2020-05-27|System and method for predicting ground layer information, and a recording medium having computer readable program for executing the method

---

CN110073246B|2021-06-11|Improved method relating to quality control

---

US9804281B2|2017-10-31|Migration velocity analysis method for VSP data

---

CN107818217B|2021-02-26|Seismic guiding method based on horizontal well exploration type visual analysis

---

Yilmaz et al.2017|Application of artificial neural networks to height transformation

---

CN110927796A|2020-03-27|Method for improving time-depth conversion precision of seismic data

---

CN107862100A|2018-03-30|A kind of shot point method for correcting error based on particle cluster algorithm

---

CN108897035A|2018-11-27|A kind of microseism weighting localization method based on wave detector weight

同族专利:

---

公开号 | 公开日

---

CN111784832A|2020-10-16|

引用文献:

---

公开号 | 申请日 | 公开日 | 申请人 | 专利标题

---

---

法律状态:

优先权:

---

申请号 | 申请日 | 专利标题

---

CN202010576436.0A|CN111784832A|2020-06-22|2020-06-22|Stratum dynamic correction method of geological model|

---